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Monitoring Wheat Powdery Mildew (Blumeria graminis f. 
sp. tritici) Using Multisource and Multitemporal Satellite 
Images and Support Vector Machine Classifier 
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Abstract: Since powdery mildew (Blumeria graminis f. sp. tritici) mainly infects the foliar of wheat, satellite remote 
sensing technology can be used to monitor and assess it on a large scale. In this study, multisource and multitempo- 
ral satellite images were used to monitor the disease and improve the classification accuracy. Specifically, four Land- 
sat-8 thermal infrared sensor (TIRS) and twenty MODerate-resolution imaging spectroradiometer (MODIS) temper- 
ature product (MOD11A1) were used to retrieve the land surface temperature (LST), and four Chinese Gaofen-1 
(GF-1) wide field of view (WFV) images was used to identify the wheat-growing areas and calculate the vegetation 
indices (VIs). ReliefF algorithm was first used to optimally select the vegetation index (VIs) sensitive to wheat pow- 
dery mildew, spatial-temporal fusion between Landsat-8 LST and MOD11A1 data was performed using the spatial 
and temporal adaptive reflectance fusion model (STARFM). The Z-score standardization method was then used to 
unify the VIs and LST data. Four monitoring models were then constructed through a single Landsat-8 LST, multi- 
temporal Landsat-8 LSTs (SLST), cumulative MODIS LST (MLST) and the combination of cumulative Landsat-8 
and MODIS LST (SMLST) using the Support Vector Machine (SVM) classifier, that were LST-SVM, SLST-SVM, 
MLST-SVM and SMLST-SVM. Four assessment indicators including user accuracy, producer accuracy, overall ac- 
curacy and Kappa coefficient were used to compare the four models. The results showed that, the proposed SMLST- 
SVM obtained the best identification accuracies. The overall accuracy and Kappa coefficient of the SMLST-SVM 
model had the highest values of 81.2% and 0.67, respectively, while they were respectively 76.8% and 0.59 for the 
SLST-SVM model. Consequently, multisource and multitemporal LSTs can considerably improve the differentiation 
accuracies of wheat powdery mildew. 
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1 Introduction 


Powdery mildew (Blumeria graminis f. sp. trit- 
ici) can occur at all stages of the wheat growth. It is 
a serious disease in some provinces of China, such 
as Sichuan, Guizhou, and Yunnan". In recent years, 
it has become more severe in the wheat-growing ar- 
eas of Northeastern China, North China, and North- 
western China. When infected by the disease, some 
serious results will be caused such as early wither- 
ing of leaves, decrease of panicle number, and de- 
crease of 1000-grain weight. Generally, the disease 
can lead to a 5%— 10% reduction in yield and when 
it is seriously infected, a more serious loss of more 
than 20% can occur”. In view of the harmful effects 
on wheat production caused by the wheat powdery 
mildew, it is of great significance to improve the 
monitoring efficiency and accuracy. However, it is 
difficult for traditional in-situ sampling and random 
investigation methods to meet the needs of large- 
scale monitoring due to the limitations in terms of 
timeliness, economy, and accuracy™. Fortunately, 
the modern information technology has facilitated 
the accurate and efficient identification of crop dis- 
eases to ensure food security” “. 

In recent years, the development of remote 
sensing technology has provided an important 
means for monitoring and forecasting large-scale 
wheat diseases and insect pests. It obtains crop in- 
formation quickly, accurately and objectively. Many 
scholars have studied remote sensing-based moni- 
toring of wheat diseases and insect pests by using 
ground-based, airborne and spaceborne remote sens- 
ing data. For example, Huang et al. showed that 
the photochemical reflectance index (PRI) was 
strongly correlated with wheat yellow rust and the 
coefficient of determination (R’) could reach 0.97 
for the PRI-based monitoring model. Zhang et al.” 


built a discriminant model for wheat powdery mil- 


dew severities by introducing continuous wavelet 
analysis based on the leaf-scale hyperspectral data. 


Luo et al.” 


constructed a monitoring model of 
wheat aphid in two-dimensional feature space de- 
rived from the modified normalized difference wa- 
ter index (MNDWI) and land surface temperature 
(LST) based on Landsat-5 TM imagery, which had 
a high discriminant precision. Zheng et al. "” con- 
structed a red edge disease stress index (REDSI) us- 
ing three red-edge bands from Sentinel-2 satellite 
imagery for monitoring the stripe rust on winter 
wheat and got a satisfying result. The above studies 
showed that the spaceborne remote sensing images 
have greatly facilitated the monitoring and diagno- 
sis of wheat diseases. 

In previous studies, multispectral satellite im- 
ages were adopted to investigate the large-scale dis- 
ease occurrence. Since the rapid development of hy- 
perspectral remote sensing technology, some studies 
on wheat powdery mildew had been explored using 
the hyperspectral data. More vegetation indices 
were derived from the hundreds of spectral bands of 
hyperspectral data. He et al."" improved the moni- 
toring accuracy of wheat powdery mildew severity 
by selecting suitable observation angles and devel- 
oping a novel vegetation index (VI). Zhao et al." 
identified the powdery mildew severities of wheat 
leaves quantitatively on hyperspectral images and 
image segmentation techniques. Khan et al."’ de- 
tected the disease using a partial least-squares linear 
discrimination analysis and the combined optimal 
features (i.e., normalized difference texture indices 
(NDTIs) and VIs). 

It can be found that the disease monitoring 
mainly depends on the vegetation changes between 
healthy and diseased wheat. Nevertheless, the com- 
mission and omission errors are usually caused due 


to the influences of other stress types such as 
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drought, inadequate nutrition, and other diseases. 
The incidence of wheat powdery mildew is in- 
volved in several affecting factors such as tempera- 
ture, humidity and planting system. Considering the 
availability of temperature data of satellite images 
for large-scale monitoring, in this study, particular 
attention was given to the contribution of LST to 
the disease occurrence. As a key habitat factor, LST 
was included in the construction of monitoring mod- 
el for the disease. LST remote sensing image was 
usually used in the previous monitoring of the dis- 
ease"), but it has a cumulative effect on wheat pow- 
dery mildew. A single-phase LST image cannot ac- 
curately represent the disease occurrence condition 
during the whole growth period. The primary objec- 
tive of this study was to explore the availability and 
feasibility to identify wheat powdery mildew using 
a combination of multisource and multitemporal 
spaceborne remote sensing imagery. More specifi- 
cally, three types of satellite images were adopted to 
identify the wheat-planting areas and retrieve the 
LST. Single and multitemporal LSTs were input in- 
to the support vector machine (SVM) to compare 


the monitoring effects. 


2 Materials and methods 


2.1 Study area 


The study area is located in Jinzhou City, He- 
bei province, China (114.97° —115.20° E, 37.80° — 
38.17° N) (Fig. 1). It has a warm-temperate conti- 
nental monsoon climate, with a flat and open ter- 
rain. There is a significant seasonal variation of sun 
radiation. Wheat is one of the important grain crops 
and widely planted in this area. The critical periods 
of wheat growth range from April to May. Due to 
its flat terrain, appropriate climate conditions and 


relatively single planting structure, the region is 


suitable for studying wheat powdery mildew using 
remote sensing technology''!. The historic statistical 
data also show that the occurrence frequency of 
wheat powdery mildew is high and hazard degree is 


serious!” 


. The wheat-planting areas were extracted 
by combining the elevation data and reflectance of 
the near-infrared (NIR) band of Chinese Gaofen-1 
(GF-1) by using decision tree classification tech- 
niques (Fig. 1(b)). The extraction results were com- 
pared with the statistical data of Shijiazhuang City, 
which could fulfill the accuracy requirement of re- 


mote sensing-based crop extraction. 
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Note: This map is made based on the standard map with the approv- 
al number of Ji S{[ 2020] 030 ( #% S[2020]030) downloaded from 
the standard map service website of the Ministry of Natural Resourc- 
es of the People's Republic of China, and the boundaries of the base 


map has not been modified 
Fig. 1 Spatial distribution of wheat-growing areas of the 


study area 


2.2 Data collection and pre-processing 


The data used in this study mainly were space- 
borne remote sensing images and field survey data 


of wheat powdery mildew. In-situ experiments were 
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carried out on 27 and 28 May 2014 in the central ar- 
eas of Jinzhou City. Several typical experimental re- 
gions were selected to collect the ground-truth data, 
where wheat was widely planted and wheat pow- 
dery mildew occurs frequently. A total of 69 valid 
data were acquired using the method presented by 
Yuan et al."”, The wheat powdery mildew severities 
were specified according to the Rules for the in- 
vestigation and forecast of wheat powdery mildew 
[B. graminis (DC.) Speer] (NY/T 613-2002). First- 
ly, the average severity D of diseased leaves for 
the colony leaves was calculated according to 
Equation (1). Then, the disease index / could be de- 
rived from Equation (2). The five levels were final- 
ly obtained according to Table 1. To increase the 
comparability of remote sensing-based disease mon- 
itoring, the five levels of severities were further di- 
vided into three levels of 0 (healthy), 1 (mild) and 2 
(severe). 


Table 1 Rules to determine the severity levels 


Level 1 2 3 4 5 
I <10 10< 7 <20 20< 7 <30 30<7<40 >40 


p- 24% 1) x 199 (1) 


where d, is the value of eight severity levels, which 
is 1%, 5%, 10%, 20%, 40%, 60%, 80%, 100%; Z, is 
the number of diseased leaves for each level; L is 
the total number of investigated diseased leaves. 

I=F XD X 100 (2) 
where F is the percentage of diseased leaf, which is 
ratio between the diseased leaves and total investi- 
gated leaves. 

Multisource remote sensing data were com- 
posed of GF-1, Landsat-8 and MODerate-resolution 
imaging spectroradio-meter (MODIS) data products. 

GF-1: Four scenes of GF-1 wide field of view 
(WFV) images were collected on 6, 17 and 26 May, 
and 7 June 2014, respectively, with the path/row 


numbers of 3/93 and 4/92. There are four spectral 
bands for the GF-1 WFV sensor within the spectral 
range of 0.45—0.89 um, with the spatial resolution 
of 16 m. The primary preprocessing procedures were 
performed including orthorectification, radiation 
calibration, atmospheric correction and image sub- 
setting. The radiation brightness (Z,, W/(m’:sr- um)) 


of a WFV image can be derived from Equation (3). 
_ DN - offset 
gain 


L (3) 


e 


where DN is the digital number of pixel. The abso- 
lute radiometric calibration coefficients (gain and 
offset) were derived from the China Center for Re- 
sources Satellite Data and Application (http://www. 
cresda. com/CN/). After finishing the radiometric 
calibration, the Fast Line-of-Sight Atmospheric 
Analysis of Spectral Hypercubes (FLAASH) mod- 
ule in ENVI 5.3 software was used to complete the 
atmospheric correction. The geometric correction 
was carried out using the second-order polynomial 
model with the accuracy of less than a pixel. 

Landsat-8: The Landsat-8 operational land im- 
ager (OLI) data were used to estimate the land re- 
flectance and the thermal infrared sensor (TIRS) da- 
ta were used for retrieving the LST"®. Four scenes 
of cloud-free images were selected and their acquisi- 
tion date was 4 April, 23 April, 15 May and 22 May 
2014, respectively. They were preprocessed includ- 
ing accurate geometric correction and radiation cor- 
rection. 

MOD114A1 product: It provided daily per-pix- 
el land surface temperature and emissivity 
(LST&E) with 1 km spatial resolution”. A total of 
20 images were obtained from the land processes 
distributed active archive center (LPDAAC) from 1 
April to 27 May, 2014. The data were preprocessed 
using the MRT (MODIS Reprojection Tool). 
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2.3 Selection of vegetation indices wheat powdery mildew were selected in this study 


i ; a based on the previous studies (Table 2). 
A total of eight VIs suitable for monitoring the 


Table 2 Selected vegetation indices in this study 


Vegetation index Abbreviation Formula 
Simple Ratio Index!” SRI SRIF Ryr/Rrea (4) 
Normalized Difference Vegetation Index 7°! NDVI NDVI= ( Ryr ~ Rrea Ra T Rea) (5) 
Normalized Difference Greenness Index”! NDGI NDGI = ( Ry 7 Rereen |N Rur kj Reveal (6) 
Soil-Adjusted Vegetation Index’! SAVI SAVI =1.5( Rug $ Racal) | Bu + Rreg + 0.5) (7) 
Enhanced Vegetation Index!” EVI EVI =2.5( Ru z Rasa} Ry + 6Rgea — 7.5Rpue + 1) (8) 
Triangular Vegetation Index™™™ TVI TVI = 0.5 120( Rye - Resa) - 200( Rrea Z Rsa] (9) 
Differential Vegetation Index'*"! DVI DVI=Ryr — Rpeg (10) 
Structure Insensitive Pigment Index'*?! SIPI SIPI = (Rii -= Rene JA Row Roi) (11) 
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Note: Rar, Ry... Ro and R 


wir? Rea Roreenand Reine 

To find out the most sensitive features, the Re- 
liefF algorithm was used due to its advantages of 
dealing with multi-class classification problem and 
having no restrictions on the data types. A sample R 
was randomly taken from the training sample set 
each time, and then k nearest neighbor samples 
(near Hits) of R were found from the sample set of 
the same class as R, and k nearest neighbor samples 
(near Misses) were found from each sample set of 
different classes of R. The input features were 
ranked according to the weights from large to small. 
Then the correlation analysis was carried out for 
each feature and the combination with the smallest 
correlation coefficient was selected as the best com- 
bination for model construction. The superiority 
and efficiency have been illustrated in remote sens- 
ing-based classification and object recognition”’””. 
Consequently, it was adopted to perform the feature 
selection, which gave different weights to the fea- 
tures in terms of the correlations between features 
and various disease samples!. Specifically, accord- 
ing to the ReliefF algorithm, all the VI variables 
were sorted in descending order of weight, and 
eight VIs were selected with the weight of 0.075 as 


the threshold value. Then, the correlation analysis 


respectively refers to the spectral reflectance of near-infrared, red, green and blue bands of a GF-1 WFV image 


among the selected features were conducted. When 
the correlation coefficient (r) of the feature owing 
the highest weight that was greater than 0.9, it was 
eliminated, and then that of the second highest 
weight with a high r was eliminated, and so on. In 
addition, there was a close relationship between the 
disease incidence and meteorological factors such 
as temperature, precipitation, humidity, etc. The 
changes of VIs calculated at different growth stages 
also affected the sensitivity to the disease. Consider- 
ing the temporal features and accumulative effect of 
temperature, three VIs were finally selected, namely 
the SAVI on 26 May 2014 and the SIPI and EVI on 
17 May 2014. 


2.4 Estimation of LST 


Two-channel nonlinear splitter algorithm was 
adopted to increase the information amount and re- 
duce the influence of errors. The equation is shown 
as (12). 


l-eé poy 


LST = by + (b, + b, —— + by) “5 I 
l-e Ae, T-T, 
(b, H b; é bs E ) 2 =d 

b (T, = TY (12) 


where ¢ and Ae respectively represents the mean 
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and difference values of the two channels’ emissivi- 
ties depending on the land type and coverage; T, and 
T, are the observed brightness temperatures of the 
two channels; b,(i = 0, 1, 2, =+, 7) is the simulation 
dataset of various coefficients which can be derived 
from laboratory data, atmospheric parameters and 
the atmospheric radiation transmission equation. In 
order to modify the calculation precision, the coeffi- 
cient depends on the atmospheric water vapour col- 
umn™!, Fig. 2 is the LST of 22 May 2014 retrieved 
from Landsat-8 TIRS image. In general, the LSTs 
of crop planting areas are lower than other regions. 
To calculate the land surface emissivity, the weight- 
ing method of vegetation coverage was adopted, 
which bases on the NDVI and vegetation coverage 
retrieved from the visible and NIR bands of Land- 
sat-8 imagery '"!, 
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Fig. 2 The LST of Jinzhou estimated by the Landst-8 TIRS 
image on 22 May 2014 


The occurrence and development of wheat 
powdery mildew is a relatively long process during 


which the LST has a significant cumulative effect. 


For example, temperature conditions which are fa- 
vorable for wheat powdery mildew infection in ear- 
ly April will aggravate the disease. As a result, the 
disease at the end of May has a high correlation 
with previous LST. The Landsat-8 LST data at four 
growth stages were selected including the standing 
stage, jointing stage, flowering stage and milking 
stage, for the retrieval of LST and their summation 
SLST is shown in Equation (13). 


4 
SLST = Ý LST, = 20 (13) 
i=1 


where SLST denotes the cumulative effect of multi- 
temporal Landsat-8 LST data; LST is calculated 
from a single Landsat-8 TIRS image; i (1, 2, 3, 4) 
represents the four stages, and the 20 means the up- 
per limit of temperature suitable for the incidence of 
wheat powdery mildew. Similarly, the summation 
of 20 MODI1A1 (MLST) can be also obtained ac- 
cording to the Equation (14). 


MLST = SLST, - 20 (14) 
2.5 Spatial-temporal fusion of MOD11A1 
and Landsat-8 LST 


Since the four Landsat-8 LST data cannot fully 
reflect the variation trend of LST during the whole 
wheat growth, the daily MOD11A1 was introduced. 
In order to obtain a data sequence that have enough 
spatial and temporal resolutions, the widely applied 
spatial and temporal adaptive reflectance fusion 
model (STARFM)"*! was used to conduct the spa- 
tial-temporal fusion of Landsat-8 LST and 
MODIIA1. The algorithm ignores the spatial posi- 
tion registration and atmospheric correction errors, 
so the pixel values of low spatial resolution (LSR) 
remote sensing data at the moment ¢ can be calculat- 
ed using the weighted sum of that of high spatial 
resolution (HSR) data (Equation (15)). 

C= TF) X A) (15) 
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where F} is the pixel value of HSR data for the posi- 
tion i at the time t; A; shows the weight of coverage 
area for each pixel; C, represents the pixel value of 
LSR data at the corresponding time. 

The STARFM algorithm first obtains the 
MODIIAI and Landsat-8 LST data and their devia- 
tion at the time ¢. The deviation is caused by the 
systematic errors and land cover changes. Mean- 
while, the Landsat-8 LST of the time ¢, in accor- 
dance with MOD11A1 were predicted. The devia- 
tion remains constant as time changes was assumed, 
and the pixel values of Landsat-8 LST at the time t, 
are Equation (16). 

HSR (Xp Yia to) = LSR (x; Yoto) + 

HSR (xa Yote) ~ LSR (XY t) (16) 

Considering the edge effect of pixels, a cloud- 
free pixel was selected in the moving window 
which was similar to the spectrum of central pixel 
when the pixel value was calculated. The calculat- 
ing the central pixel value is shown as Equa- 
tion (17). 

HSR (Xy25 Yans to) = 

SSS Wi X (LSR (x; Yi to) + 


EENES 

HSR (x; Yi t) — LSR (x; Yi ty) ) (17) 
where w is the size of the moving window; 
(XwosVw2) represents the position of central pixel; 
W,,, denotes the weight coefficient of a pixel similar 
to the central pixel. The spectral distance weight, 
temporal distance weight, and spatial distance 
weight of a similar pixel were obtained in the win- 
dow by a normalization method. The three weight 
coefficients were taken by referring to the study of 
Gao et al. in 2006". Four Landsat-8 LST (LST) 
and 20 MOD11A1 (LST,) were used for the spatial- 
temporal fusion. The fusion data sequences were 


summed up as the SMLST in Equation (18). 


4 20 
SMLST = LST, + YLST, (18) 
i=l j=l 


2.6 Construction of monitoring models 


Support vector machine (SVM) has the advan- 
tages of a simple structure, strong generalization 
ability, and high accuracy, which has been widely 
used in the classification of remote sensing imag- 


[36] 


es”. The discriminant function of the model is 


shown as Equation (19). 

£0) =n] Sark o (19) 

“eS, 

where a, is the Lagrange multiplier; S, is the support 
vector; x; and y, represent two kinds of support vec- 
tors; b is the threshold; k(x, x) represents a positive 
definite kernel function which satisfies the Mercer 
theorem. 

The Z-score method was used to standardize 
VIs and temperature data according to Equation (20) 
due to their different units. The training and valida- 
tion datasets were divided using the ratio of 7:3. 
Four SVM-based models (LST-SVM, SLST-SVM, 
MLST-SVM and SMLST-SVM) were trained using 
the optimally selected three VIs. In addition to the 
VIs, these models were also involved in the Land- 
sat-8 LST on 22 May (LST), four cumulative Land- 
sat-8 LST data (SLST), 20 cumulative MODIIAI 
(MLST) and cumulative spatial-temporal fusion 
LST combing Landsat-8 LST and MODIIAI 
(SMLST), respectively. These models will explore a 
better application of remote sensing-based LST to 
the monitoring of wheat powdery mildew. 

x'=(x-p)/p (20) 
where x’ is the standardized data; u is the mean of 


original data; p is the standard deviation. 


3 Results and discussion 


3.1 Validation of the monitoring results 


The cross-validation was adopted to estimate 


the monitoring accuracies. As shown in Table 3, 
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four indicators, namely the user accuracy (UA), pro- 
ducer accuracy (PA), overall accuracy (OA) and 
Kappa coefficient (k) were used to assess the four 
SVM-based models. It could be seen that the SLST- 
SVM and SMLST-SVM achieved better classifica- 
tion performance. In terms of OA, the SMLST- 
SVM model obtained the best result, followed by 
the SLST-SVM and MLST-SVM models, while the 
LST-SVM model got the lowest value. The OAs of 
models indicated that the LST had a cumulative ef- 
fect on the wheat powdery mildew infection. The k 
values of SMLST-SVM, MLST-SVM SLST-SVM 
and LST-SVM were 0.67, 0.54, 0.59 and 0.38, re- 


spectively, which also showed the similar trend with 


OAs. From the perspective of UA, the ability of the 
four models to distinguish diseased and healthy 
wheat was strong. The differentiation ability of 
SLST-SVM and SMLST-SVM models for mild and 
severe wheat powdery mildew were significantly 
higher than that of LST-SVM model. The accura- 
cies of MLST-SVM model were slightly lower than 
those of SLST-SVM and SMLST-SVM models, 
mainly due to the low spatial resolution of 
MOD11A1 data on the city scale. The above results 
show that the introduction of multi-temporal and cu- 
mulative LST can effectively improve the monitor- 
ing and identification of wheat powdery mildew se- 


verities. 


Table 3 Comparison of the accuracies using the four monitoring models 


Confusion matrix 


Accuracy assessment indicator 


— Healthy Mild Severe User accuracy/% Producer accuracy/% Overall accuracy/% k 

Healthy 14 T 0 66.7 63.6 

LST-SVM Mild 8 25 4 83.8 65.8 63.8 0.38 
Severe 0 6 5 45.5 55.6 
Healthy 14 7 0 66.7 82.4 

SLST-SVM Mild 3 32 2 86.5 74.4 76.8 0.59 
Severe 0 4 7 63.6 77.8 
Healthy 15 6 0 71.4 88.2 

MLST-SVM Mild 2 32 3 86.5 TA T39 0.54 
Severe 0 7 4 36.4 57.1 
Healthy 16 5 0 76.2 84.2 

SMLST-SVM Mild 3 33 1 89.2 78.6 81.2 0.67 
Severe 0 4 T 63.6 87.5 


Huang et al."” identified wheat powdery mil- 
dew of the study area using 30 m-resolution Chi- 
nese HJ-1A/1B data to inverse LST and extract four- 
band reflectance data and build seven VIs. A combi- 
nation method (GaborSVM) of SVM and Gabor 
wavelet features were proposed to obtain the OA of 
86.7% that was higher than 81.2% of this study, but 
the OA of SVM-based method was 80% and lower 
than this study. The primary reason is that the spa- 
tial resolution of HJ-1B IRS is 300 m, but it is just 
1000 m for MOD11A1. The comparison study 


shows that spatial resolution of optical and thermal 
infrared satellite remote sensing images is an impor- 
tant factor of affect the accuracy of wheat powdery 


mildew. 
3.2 Mapping of wheat powdery mildew 


Based on multisource and multitemporal Land- 
sat-8, GF-1 and MODIS data, three VIs of SIPI, SA- 
VI and EVI were optimally selected by the ReliefF 
algorithm and correlation analysis. The VIs and 
four temperature data of LST, SLST, MLST and 
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SMLST were respectively used to construct four 
monitoring models through the SVM, namely the 
LST-SVM, SLST-SVM, MLST-SVM and SMLST- 
SVM. For example, the severity distribution on 26 
May 2014 in Jinzhou was shown in Fig. 3 using the 
SLST-SVM and SMLST-SVM models. The overall 


spatial distribution of wheat powdery mildew using 
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the two monitoring models were similar. Neverthe- 
less, there were also some obvious differences as 
shown in the red boxes. It was more serious in the 
eastern part than in the western part of the study ar- 
ea. It was also obvious that the wheat powdery mil- 
dew mainly occurred in the areas where wheat was 


widely planted. 
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Fig. 3 Comparison of monitoring wheat powdery mildew by the SLST-SVM model and SMLST-SVM model 


In the central regions where the ground-truth 
points were located, the monitoring results of the 
two models were also similar by visual observation. 
In the two figures, the distribution trends of wheat 
powdery mildew were both relatively concentrated. 
In comparison with the 32% of healthy samples, 
55% of mild samples and 13% of severe samples 
for the in-situ investigation, they were 37%, 49% 
and 14% for the SLST-SVM model and 31%, 55% 
and 14% for SMLST-SVM model, respectively. It 
can be seen that the SMLST-SVM model had a bet- 


ter result than SLST-SVM model. 
3.3 Analysis of influence factors 


Temperature is one of the key factors to affect 
the incidence of wheat powdery mildew, however, it 
is difficult and inaccurate to identify the disease in 
small and medium-sized regions, due to constraint 
of low spatial resolution of regular meteorological 
data. Ma et al.°*! combined meteorological and re- 
mote sensing data to monitor wheat powdery mil- 
dew, because the distribution of meteorological sta- 


tions was too sparse. The primary objective of this 
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study is to compare the relationship between LST 
and the disease severities of wheat powdery mil- 
dew. It can be found that different LST data have 
significant impact on the accuracies of the disease 
diagnosis. For example, in the retrieval of LST from 
Landsat-8 data, the two-channel nonlinear splitter 
used in this paper was more reliable than the single- 
channel method"’!. The spatial resolution of Land- 
sat-8 TIRS data and the temporal resolution of 
MODIIAI1 data can both meet the requirements. 
Therefore, the combination and fusion of both the 
images were performed to acquire a better perfor- 
mance. The accumulated temperature and effective 
accumulated temperature are key factors to affect 
wheat powdery mildew’. Multitemporal LST is 
better to show the temperature influence on the dis- 


ease monitoring than a single one. 
4 Conclusion 


It is a progressive process for the occurrence of 
wheat powdery mildew, the cumulative effect must 
be considered with the disease development. To ac- 
curately identify the disease, multisource and multi- 
temporal GF-1, Landsat-8 and MODI1A1 were 
jointly utilized at a city scale. The vegetation fea- 
tures and LSTs are simultaneously adopted to mon- 
itor the disease for improving the classification ac- 
curacy. 

To find out the most sensitive vegetation fea- 
tures to wheat powdery mildew, four most impor- 
tant growth stages were selected including the stand- 
ing stage, jointing stage, flowering stage and milk- 
ing stage, which are also the criteria for selecting 
satellite images. Additionally, it has a progressive 
process for the disease incidence with the growth of 
wheat and changes of meteorological factors. As 
one of the most essential influencing factors, the ac- 
cumulative effects in LSTs must be considered 


when identifying the disease. Landsat-8 TIRS has 


high spatial resolution but low temporal resolution, 
however, it is quite contrary to MOD11A1. As a re- 
sult, the STARFM was selected to perform the spa- 
tial-temporal fusion of both data. Four SVM models 
were respectively constructed through a single 
Landsat-8 LST (LST), multitemporal Landsat-8 
LSTs (SLST), cumulative MODIS LST (MLST) 
and the combination of cumulative Landsat-8 and 
MODIS LST (SMLST). 

The OA of LST-SVM model was improved by 
13% using the four temporal LSTs. In addition, the 
Kappa coefficient also increased from 0.38 to 0.59, 
indicating that LST is a key habitat factor for the oc- 
currence and development of wheat powdery mil- 
dew, due to a cumulative effect. On the other hand, 
the accuracies of MLST-SVM model were smaller 
than that of SLST-SVM and SMLST-SVM models, 
indicating that it is inappropriate to apply MOD11A1 
data directly to the monitoring of wheat powdery 
mildew in a relatively small area. Conversely, the 
performance of SMLST-SVM is slightly better than 
that of SLST-SVM, indicating that the monitoring 
performance can be enhanced by using a combina- 
tion and fusion of high-resolution spatial-temporal 


remote sensing data. 
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